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Abstract 

Inspired by the rapid elongation of air columns injected into vortices by dolphins, we 
present an exact inviscid solution for the axial speed (assumed steady) of propagation 
of the tip of a semi-infinite cylindrical bubble along the axis of a cylindrical vortex. The 
bubble is assumed to be held at constant pressure by being connected to a reservoir, the 
lungs of the dolphin, say. For a given bubble pressure, there is a modest critical rotation 
rate above which steadily propagating bubbles exist. For a bubble at ambient pressure, 
the propagation speed of the bubble (relative to axial velocity within the vortex) varies 
between 0.5 and 0.(i of the maximum rotational speed of the vortex. Surprisingly, the 
bubble tip can propagate (almost as rapidly) even when the pressure minimum in the 
vortex core is greater than the bubble pressure; in this case, solutions exhibit a dimple 
on the nose of the bubble. A situation important for incipient vortex cavitation, and 
one which dolphins also demonstrate, is elongation of a free bubble, i.e., one whose 
internal pressure may vary. Under the assumption that the acceleration term is small 
(checked a posteriori), the steady solution is applied at each instant during the elonga- 
tion. Three types of behavior are then possible depending on physical parameters and 
initial conditions: (A) Unabated elongation with slowly increasing bubble pressure, and 
nearly constant volume. Volume begins to decrease in the late stages. (Bl) Elongation 
with decreasing bubble pressure. A limit point of the steady solution is encountered 
at a finite bubble length. (B2) Unabated elongation with decreasing bubble pressure 
and indefinite creation of volume. This is made possible by the existence of propagat- 
ing solutions at bubble pressures below the minimum vortex pressure. As the bubble 
stretches, its radius initially decreases but then becomes constant; this is also observed 
in experiments on incipient vortex cavitation. 


1 Introduction 

Marten et al. [1] report captive dolphins making bubbles interact with vortices in a number 
of ways. This behavior is performed, apparently for play, without any training or food 
reward, and at the dolphins’ own volition. 

In one case, the dolphins make ring vortices with a flick of their tail fins and then 
turn around to inject air into them from their blowholes. The air forms a tubular bubble 
that lengthens rapidly along the vortex core, eventually forming a toroidal bubble; see 
photographs at the top of pp. 86-87 in Marten et al. [1]. Dolphins also inject air into the 
longitudinal vortices created behind their dorsal fin. This too forms a cylindrical air column 
that rapidly propagates along the vortex; see figure 1. This situation differs from the case 
of the vortex ring in that axial flow is present in the vortex. In both of these examples, the 
bubble remains connected to the dolphin’s blowhole or mouth while undergoing substantial 
increase in length. Hence we shall initially assume that the bubble pressure equals the 
reservoir pressure, which in turn is taken to be constant in time. 

The mechanism for propagation of the bubble is simply that the lower than ambient 
pressure in the vortex core sucks the air, provided the surface tension is not too high. To 
see this mathematically, consider figure 2 which shows one end of the bubble propagating 
into the vortex. The bubble is assumed to be sufficiently long that there exists a region 
(station 2 in figure 2) where axial derivatives may be neglected; this is made more precise in 
§4. In order to the make the flow steady between stations 1 and 2, choose a reference frame 
that moves with the speed U b of the bubble nose. Constancy of Bernoulli head along the 
axis streamline implies that the pressure at the tip of the nose is p n = p a + p(U ~ Ub) / 2 > 
where p a is the pressure at the axis of the undisturbed vortex, p is the liquid density, and 
U is the axial velocity of the liquid at the vortex axis in the laboratory frame. Combining 
this with the pressure matching condition p b - Pn = gives the bubble speed: 

U h = U± ^(Pb - Pa - Tk), (!) 

where T is the surface tension and k is the total curvature at the tip. The ± symbol 
means that the tip may either advance through the vortex or retract. We shall choose the 
physically observed advancing solution; in reality it is presumably the initial condition that 
sets up the elongation. The speed would be determined by (1) if only k were known. In 
the present work, one is fortunate to be able to evaluate the speed without knowing k, and 
in fact, to able to deduce n a posteriori using (1). What makes this evaluation possible is a 
useful method of analysis introduced by Batchelor [2] in which one matches two cylindrical 
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(axially invariant) solutions of the steady Euler equations at stations 1 and 2. 

The case of a free bubble not connected to a reservoir and therefore having variable 
internal pressure is also important in three contexts: 

(i) First, the dolphin stops blowing at some point. 

(ii) Dolphins have also been observed to release free spherical bubbles and then to “whip 
them into shape” using vortices. For instance, one dolphin released a cloud of bubbles from 
its blowhole while swimming rapidly; the bubbles elongated into a long tube when they 
encountered the longitudinal vortex behind the dorsal fin; see photographs at the bottom 
of pp 86-87 in Marten et al. [1]. One observes a tremendous creation of volume from a few 
small bubbles. This may seem contrary to intuition for if the volume increases, the bubble 
pressure must at some mstant fall below the vortex pressure, thus eliminating the driving 
mechanism for the elongation according to (1). The present solutions indicate that what 
permits elongation in this situation is k being < 0. Another dolphin stood upside down m 
the water and released a bubble from its mouth which rose towards the surface; when the 
bubble reached the tail fin, it was elongated into a ring vortex by a flick of the tail fin. 

(iii) A similar phenomenon was observed in experiments on incipient vortex cavitation 
behind a wing tip conducted by Maines and Arndt [3]. They observed four stages: (I) 
Motion of a small spherical bubble towards the vortex axis driven by the radial pressure 
gradient in the vortex. This is also the reason why the bubble is initially helical in figure 1. 
(II): Volume expansion of the bubble in the shape of an ellipsoid whose major and minor 
axis both grow, with the major axis growing much more rapidly. (Ill) Explosive growth 
of the major axis accompanied by a decrease in minor axis. (IV) In the final phase, the 
bubble grows as a cylindrical bubble in which the length continues to increase while the 
radius becomes constant. 

We shall provide a treatment in which the steady solution is applied to such situations 
under the assumption of quasi-steadiness. 


2 Analysis 


2.1 Bubble with specified pressure 

We employ cylindrical coordinates ( x , cr, <j>) with corresponding velocity components (u, v, w). 
The streamfunction ip(x,cr) defined by 


u — 


1 d'lp 
a da ’ 


v = — 


1 d'lp 
a dx ’ 


( 2 ) 
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Figure 1: Dolphin injecting air into a vortex it previously created using its dorsal fin. 
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Figure 2: Sketch used for the analysis. 
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describes the flow in a meridional (xa) plane. The swirling flow in the undisturbed vortex 
at station 1 is taken to be a Rankine vortex 

_ f flcr, a < 5i; (3) 

1 | fi(5 2 /cr, a > <5i, 

where Si is the vortex core size. The axial velocity in the undisturbed vortex is assumed to 
have a uniform axial velocity, U, such that in the reference frame moving with the bubble 

tip we have 

_ f U-U b , a <6i] (4) 

\ ~U b , a > Sl 

These choices for the vortex flow were made for analytical tractability. More general profiles 
may be treated by the procedure below but it has to be implemented entirely using numerical 

analysis. 

Steady axisymmetrie inviscid flow is governed by the Squire-Long equation [2], 

1 d± = ff 2 H'{tp)/p - K{iP)K'{iP), (5) 

da 2 da 2 a da 

where K(ip) = aw is the circulation function and H{i>) = p + p|u| 2 /2 is the Bernoulli head. 
We shall follow a very useful but simple method of analysis introduced by Batchelor [2] m 
which one joins two axially invariant solutions of (5) at two stations. 

The first step is to obtain the functions K(ip) and H{ip) from the flow at station 1. At 
this location, the streamfunction obtained using (2) and (4) is. 


i/>i(cr) 


(1 U - U h )a 2 / 2, cr < *i; 

{U8\-U b a 2 )l 2, a > Si. 


For later use, let us define the streamfunction for the vortex boundary: 


■ l p 0 = t j; 1 (8 l ) = (U -U h )8l/2. 


(6) 

(7) 


After solving (6) for cr 2 in terms of if), the circulation function is obtained. 


K{ip) = aw = 


2ipnf(U-U b ), a <5 1; 
nsj, a > Si . 


The pressure at station 1, obtained after integrating 

dpi _ pw{ 
da a 


( 8 ) 

(9) 
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Poo + P n 2 (a 2 /2 - «?) , cr < ^i; (io) 

Poo~P^i/(2a 2 ), C>S V 
Evaluating Hty) = Pi + p|ui| 2 /2 and writing a in terms of i> gives 

I Poo + P& (2 */(U - u b ) -6 2 )+p(U- U b ) 2 /2, within the vortex; (u) 
W) “ \ Poo + pU b /2, outside the vortex - 

Inserting (8) and (11) into (5) gives for an axially invariant flow: 

d 2 ip 
da 2 

We obtain V at station 2 by solving (12). In the interior region, (12) has the same form 
as an equation obtained by Batchelor [2](p. 546) and we simply use his solution. In the 
exterior, we assume a solution of the form ip oc a m . The result is. 


1 dip 2ft 2 J a 2 - 2 ip/(U - U b ), inside the vortex; 

a do ~~ (U - U b ) | 0, outside. 


(12) 



_ ([/ b _ jj) a 2 ! 2 + a [AJ\(ko) + BYi(ka )} , a < 8 2 ; 
-U b a 2 /2 + C, 


(13) 


where 


k = 


2 ft 


(14) 


(t/ - f/b) ' 

The constant C = (l/2)U5'( by imposing the condition that ip 2 {8 2 ) = V’O- The velocity 
corresponding to (13) is 

- (f/b - U) + AkJo(ka) + BkY 0 (ka), a < <5 2 ; / 15 j 

-f/b, a ^ 52 ‘ 


u 2 (cr) s: 


The five unknowns U b - U, <5 b , ,A,B are obtained by imposing five conditions: 

(1) The bubble surface must be a streamtube, i.e., ip 2 (S b ) - 0: 

-\{U b -U) + AJi(k6 b ) + BYi{kS b ) = 0. (16) 

(2) The vortex boundary must be the same streamsurface at stations 1 and 2, i.e., ip 2 {8 2 ) — 

ipo- i 

-- (U b - U) 6 2 = -- (f/b -U)6 2 + S 2 [AJi{kd 2 ) + BYi{kS 2 )] . (17) 

2 ^ 

(3) The pressure must be continuous at a = 6 2 , i.e., p 2 (a = S 2 ) = p 2 {o = <^)- The P ressure > 
p 2 (a = S 2 ), just outside the vortex is obtained by equating the value of the Bernoulli head 

in the exterior region to its definition: 

Poo + ^pt/ b — P?( a = ^2) + 2 P [ w 2 (^) + ^2(^2")] • ( 1 ^) 
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The pressure, pa(ij), just inside the vortex is obtained by evaluating H{i>) of the inner 
region at ^ &nd equating it to its definition. 

Poo + ~ = ^ 2 (^ 2 ) + 2 P [ U 2^2 ) + w 2^2 )] • 

Subtracting (19) from (18) and using the fact that both pressure and azimuthal velocity 
must be continuous gives 

“2(^2 ) = ± {Ub - u) ( 2 o) 

To choose the correct sign in (20), one argues that in the case when the axial velocity in 
the vortex U = 0, a vortex sheet does not exist and therefore the axial velocity must be 
continuous, i.e. u(<5 2 -) - -U h . Hence, we choose the minus sign in (20) and obtain: 

AJ 0 (k8 2 ) + BY 0 {k8 2 ) = 0. (21) 


(4) Forces must balance on the bubble surface at station 2. 

p h -p 2 (S h ) = T/6 h . (22) 

The pressure p 2 (<5b) is obtained by equating the value at the bubble surface {i> — 0) 

to its definition: 

Poo - ptfsl + l -p{U h - U ) 2 = p 2 (S h ) + \p [«l(*b) + ti|(«b)] • (23) 

Solving (23) for p 2 (^b)> substituting this into (22), and noting that the azimuthal velocity 
w(5 b ) on the bubble surface vanishes because K{ip = 0) = 0 gives: 

Pb ~ Poo ~ 2 P (^ b - U) 2 + l + 2^(^b) ~ ^ = 


where 

ri;(5b) = ~(U\) “ U) + AkJ§(kS\ j) 4- BkY{}(k5 b)* (25) 


At this point, let us introduce non-dimensional quantities (denoted by hats) using the 
density p, vortex core radius <5i, and angular velocity as the scaling quantities. 


A 

Apb 


4- B 

Pb ~ Poo 
pfl 2 5 2 ’ 


B 

fW 2 

f = 



T 


pVt 2 5\ 



A U h = 


UbzJL 

fWi ’ 


(26) 


An expression in dimensional form may be converted to non-dimensional form by setting p, 
<5i, and Q. to unity and placing hats over all quantities. 
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(5) The fifth condition in that axial momentum must balance in the region between stations 
1 and 2: ^ 

27rT5 b + 2 tt J ado [pi(a) + pi^(er) - p 2 {o) - i( CT )j = 0 ^ 

Notice that fpuvodx, the radial flux of axial momentum at the outer radial boundary 
oo, does not appear in (27). This is because the radial velocity v ~ a 2 given that 
there are no sources of volume. When pressure and velocity profiles are substituted into 
(27), several integrals involving Bessel functions appear. These integrals are provided in 
Gradshteyn & Ryzhik [4] 5.54.2, 5.56.2, and Watson [5] p.134, eq. (10). The integrals were 
supplied to Mathematica which performed the rest of the algebra. The final result is: 

Ti+T 2 +T 3 =0, ( 28 ) 


where 


Ti 


T 2 


t 3 


3 1 -J2 1 J2 , 1 1 ?4 _ 1 

"8 + 2 <5 ' 2_ 2 Jb + 8 5b ~ 8* 2 2 

^[A 2 (Sg/i -S 2 f2) + AB (P b h 

U h (a/,. + Bfs) 


6ip h + 6 h T + ^6lu£ - l - log (j 2 ) 

-P 2 h)+B 2 (P h h-P 2 h)} 


with 

fl = JoiVb) + Joiv^hiPb) 

/ 2 = Jo(m) + Jo{V2)J2{V2) 

— J 2 (%)^o( , ?b) 4" 2Jo ( ? 7b)^o( T ?b) A Jo{' l 7b)Y 2 {T]b) 

/ 4 = 2J 0 {rj 2 )Yo{r) 2 ) + Mm)Yo{V 2 ) + MV 2 )Y 2 {m) 

/.5 = Y o{ r ib) + yo(77b)y>(%) 

h = >oM + Yo(m)Y2(m) 

fl = foJlim) - tbJliVb) 

fg - <5 2 Yi(t? 2 ) - £ b Yi(77b) 

in which 2 

T ] b = k <5 b , m^kfa, k = ( 29 ) 

It is not difficult to show that energy conservation is automatically satisfied by the 
requirement that H(f) is the same function of ip at both stations. Similarly, angular 
momentum is conserved because K {ip) is the same function everywhere. 
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The five equations (16), (17), (21), (24), and (28) were solved using routine NEQNJ in 
the IMSL library which uses a variant of Newton’s method. Elements of the Jacobian matrix 
were analytically calculated with the help of Mathematica. We initially found it impossible 
to find a starting guess for the five unknowns that would converge to any solution. To 
overcome this difficulty, we pretended that the problem was that of a solid cylinder of 
specified radius being pushed through a vortex at a specified speed. Only the first three of 
the five conditions are then needed and the method readily converged to a solution when 
starting with an initial guess analytically obtained for the small swirl limit. This converged 
solution was then substituted into the remaining two conditions (which happen to be linear 
in T and Ap b ) to obtain the surface tension T and bubble pressure Ap b that would make the 
solution valid for a bubble. Starting with this solution Eve, all the other solutions presented 
in this paper were obtained by gradually varying parameters and using the previous solution 
as a starting guess. 

Since U always occurs in the combination U h - U, the effect of a uniform axial flow on 
bubble speed is merely additive. The five equations remain unchanged when A, B and AC/ b 
all change sign. Thus, for every solution that represents a bubble advancing into the vortex, 
there exists is a solution that represents a bubble withdrawing from the vortex. This is a 
consequence of time-reversal symmetry of inviscid motion and in practice initial conditions 
determine the relevant solution. 

2.2 Free bubble 

Consider a cylindrical bubble at time t with finite length L(t) and internal pressure p h {t). 
If the gas within the bubble remains adiabatic, then 

Pb (t)Vi(t)=p b0 V 0 \ (30) 

where the subscript ‘0’ denotes initial values. The bubble volume V{t) « nSl(t)L{t) if 
the bubble is long enough that the contribution from its ends is negligible. Writing p b = 
Poo + Ap b and p b o = P<, o + Ap b0 , (30) gives 

l= /poo + Apboy 7 ^ (31) 

V Poo + Apb / 1x8^ 

If the acceleration Ou/dt is small compared to the pressure gradient, say, then one may 
make the quasi-steady assumption that at each instant, the bubble radius and speed of 
extension of each end can be obtained from the steady solution at the bubble pressure 
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implied by (31). We estimate: 


du d(A U b ) d(A U b ) dL _ OArr d{AU b ) _ d(AU& 
&t ~ ~~dT~ ~ dL dt - l * Uh dL dL ’ 


and 


VP~ 


pLt 2 8 2 

~T~‘ 


In non-dimensional terms therefore, the requirement for quasi-steadiness is: 


(33) 


_ (dAU 2 ) 


c qs — 


dL 


< 1 . 


(34) 


Let us consider two cases for the ‘O’ state: (i) a spherical bubble of radius R b in quiescent 
liquid that is then subjected to the vortex flow field, and (ii) a cylindrical bubble just after 
the injection device closes. In the first case, Apbo = 2T/f?b and Vo = (4/3 )nR^. In 
the second case Apbo is the reservoir pressure, and Vo = ir8^ 0 Lo where <5bo is obtained 
from the steady solution for the given Api,o ■ Substituting these into (31) we can calculate 
L corresponding to any solution point of the steady problem. Note that we cannot say 
anything about the path from the ‘O’ state to the cylindrical one: for instance, a dynamical 
barrier may exist that prevents the latter from being realized. 


3 Results 

3.1 Bubble with specified pressure 

Figure 3 displays various results as the relative bubble pressure is varied for different values 
of the surface tension. Focus first on the propagation speed (graph a) and notice that two 

of the curves ( and ) have a break near Apb = 0. This is where the solution 

algorithm stopped converging either because the solution does not exist or changes too 
rapidly to track. 

Each curve also has a limiting value of Ap b below which the algorithm also stops con- 
verging. We believe these to be genuine limit points beyond which steady solutions do not 
exist. The presence of a limit point is partially plausible from (1) which shows that the 
radicand becomes negative if Ap b is sufficiently large. However, (1) implies that AU b = 0 
at the limit point provided n continues to exist. The fact that the actual limit point occurs 
for higher AU b means I hat other constraints prevent a solution from existing. 

Figure 36 plots the vortex radius at station 2 and shows that the vortex pinches instead 
of being displaced outward by the bubble. 
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Figure 3: , ° : T — 1; , A : T — 0.75; , a : T — 0.40; : T 0.10; 

: f = 0.05; : f = 0.001. (a) Bubble propagation speed A U h . (b) Disturbed 

vortex radius ((f 2 )- ( c ) Bubble radius 6^. (d) Curvature index a of the bubble tip. 



Apb 


Figure 4: Bubble radius compared approximate formulas. Each curve is labelled with the 
value of T. , exact; , equation (36); , equation (37). 
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Figure 3c plots the bubble radius. One may obtain an approximate formulas by sup- 
posing that the flow at station 2 is negligibly disturbed relative to the vortex at station 1. 
The equilibrium radius of a cylindrical bubble in such a flow is given by the condition: 

Pb~Pi(S b )=T/S h , (35) 


(36) 


(37) 


which, upon using (10), becomes in terms of non-dimensional variables: 

+ 1 = T/8 b . 

An even simpler approximation is obtained by neglecting the term quadratic in <$b : 

? T 

b Ap b + 1 

Figure 4 (dotted) shows that (37) is slightly better than (36) and that both approximations 
overpredict the exact result implying that the liquid pressure at the bubble radius is smaller 
than at the same radius in the undisturbed vortex. 

Figure 3d plots the curvature of the bubble tip as deduced from (1) and the fact that 
p a = Poo - pQSi- Actually, a curvature index a, defined as 


■ = 1 — kSu* 


(M) 


is plotted in order to emphasize small difference from kS\> = 1 on a log scale. A hemispherical 
nose has k — 2/S b {a = -1), a flat tip has k = 0 (a ='l), while a nose with negative 
curvature has k < 0 (<* > 1). Figure 3d shows that not only is the tip is flatter than 
hemispherical, i.e., ct > —1 but ct > 0. That this should be so is proven as follows. 
Constancy of the Bernoulli head from the tip to a point on the bubble surface at station 2 
gives: 

Pn = P 2 (^b) + 2(^2) = Pb - T/ 5 b + ” ^2(^2), ( 3 ®) 

where the fact that the azimuthal velocity is zero on the bubble surface has been used. On 
the other hand 

Pn = Pb “ kT. (40) 

Equating (39) and (40) gives K.5 b < 1. 

A word of caution: t he fact that a « 0 for large Apb i n figure 3d should not lead one to 
say that the curvature is approximately known and therefore (1) may be used to approximate 
the speed. Substituting into (1) the fact that the bubble radius 5 b « T/(Ap b + 1) for large 
Ap b gives: 

AU b « \f2a (Apb + 1). (41) 
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Figure 5: Propagation speed of a bubble at ambient pressure ( Apt, = 0) as a function of 
surface tension. 


Thus for large Ap^ the speed involves the product of a, no matter how small, and a large 
Apb 

Note from figure 3a that when T is small enough, the bubble is able to propagate even 
when Apb < — 1, i.e., when bubble pressure is less than minimum vortex pressure. In this 
case the tip curvature is negative, i.e., the nose is dimpled (figure 3d). 

The case of Apb = 0 is relevant for an unpressurized reservoir and figure 5 plots the 
speed as the surface tension is varied. Solutions exist for T < .380 which for air and water 
and a 3 cm core radius corresponds to a modest ST/(2-7r) > 0.9 rev/s. 

3.2 Free bubble 

Here we consider the quasi-steady behavior of a free cylindrical bubble whose pressure and 
dimensions are compatible with it having initially been a spherical bubble of radius Rq in 
a quiescent liquid. The parameter space now includes Rq = Ro/&\ and p <*, = Poo/(pH 2 ^i) 
in addition to T and Apb- Instead of mapping out the space, we shall be content to sketch 
three scenarios which we shall denote as A, Bl, and B2. 

Consider (31) (placing hats on all quantities) as Apb varies. If we begin in the regime 
Apb •C Poo i then the first factor in (31) is practically constant as Apb varies. Since, the 
bubble speed is positive, L must increase and therefore 6^ must decrease. Two cases are 
then possible: the cylindrical bubble may be at a point in figure 3c where the <5b versus 
Apb curve either has a (A) negative or (B) positive slope. 

In the former case, the solution point must move to the right along the curve to allow 
<?b to decrease. Thus Apb increases and the elongation continues unabated. The bubble 
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volume decreases with elongation, however, as long as Apb ^ Poo? the decrease is negligible 
and the bubble will be nearly incompressible. At some point during the elongation, A p b will 
become comparable to poo- This means that the first factor in (31) will decrease and so 5b 
must decrease faster to allow L to increase. At this stage, the bubble volume will begin to 
wither away. In fact, when Apb ^ Poo we may use (37) to conclude that the bubble radius 
will diminish as 5 h ~ L - 2 7/( 2 7-i) - L -1 - 56 , which is faster than the 8 h ~ £~ 1/2 behavior of 
an incompressible bubble. 

For case Bl, the 8 h versus Ap b curve has a positive slope. Hence the solution point 
must move to the left along the curve: Apb becomes more negative and eventually a limit 
point is encountered which presumably stops the elongation. 

Case B2 is similar to Bl except that the location of the limit point Apb,lim < — Poo 
allows unceasing elongation. As Apb becomes more negative during the elongation, it will 
become comparable to -poo &t some point. The first factor in (31) then increases and it is 
no longer necessary for the bubble radius to shrink to allow L to increase. At this point, the 
bubble radius becomes constant and volume begins to be created indefinitely with bubble 
pressure becoming arbitrarily close to vacuum. 

We now illustrate each of the three cases by means of specific examples. 

To demonstrate type A behavior, let us initially have a spherical air bubble of 2 mm 
radius surrounded by water at atmospheric pressure, that is subjected to a vortex with a 2 
cm core rotating at 2 rev/s. These values imply T = .058, Poo — 1604, and Rq —*0.1. The 
first step is to generate a solution locus for T = .058 by varying Apb us in figure 3. We 
find that the locus has two branches similar to the dotted curves in figure 3. The present 
analysis cannot say which branch the bubble will end up on. The left branch exhibits type 
Bl elongation terminating quickly at L = 0.39. The neglect of axial derivatives is unlikely 
to hold for such a small L (see §4) ; still behavior with similar characteristics may exist in 
reality. The right branch exhibits unassailable type A behavior as shown in figure 6: the 
pressure increases with length and the elongation continues unabated. Figure 6 b shows 
that the bubble radius diminishes very nearly as L -1 / 2 , i.e., the bubble remains very nearly 
incompressible. As noted above, we should expect a stage to be reached (when Apb > 1604) 
where volume diminishes and its existence was verified. The measure of merit, e qs , of the 
quasi-steady assumption, decreases from 0.2 to 0.004 in the range of figure 6. 

To depict type B elongation, let us choose f = 0.001 (see chain-dot curves in figure 3). 
This is achieved by increasing f)/(27r) to 15.18 rev/s, leaving all other dimensional param- 
eters the same as in the previous example. We then have Poc = 27.83. Again, the solution 
locus has two branches The right-hand branch (type A) is found to start at L ~ 500 and 
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Figure 6: Evolution of bubble pressure and radius with elongation (example illustrating 
case A). 


one cannot rule out an unsteady process that somehow produces such a cylindrical bubble. 
Let us focus on the left-hand branch (type B). Figure 7 a shows that the bubble pressure 
decreases as the bubble elongates, and a limit point is suddenly encountered. The curves 
were truncated at the left in order that the acceleration parameter 2 x 10~ 5 < e qs <0.1 
in the range of the plot. Again, quasi-steadiness becomes a better approximation as the 
bubble extends. 

Finally, to obtain type C behavior, let us choose fi/(27r) = 10 rev/s, S\ = 7 cm and 
Ro = 3 mm. We then have T — 5.4 x 10” 5 , Poc = 5.24, and Rq = .043. The right branch 
(type A) begins at an absurd L « 14,000 and so we focus on the left branch which begins 
at a reasonable state. The limit point of the left branch lies at Apb,iim ~ —5.42 < —poo thus 
allowing unabated elongation together with creation of volume. Figure 8a shows that the 
bubble radius initially decreases and then becomes constant as seen in stages III and IV of 
the experiments [3]. The bubble pressure, shown in figure 8a , becomes nearly constant as 
the bubble elongates. This means that the bubble speed becomes nearly constant and the 
acceleration parameter becomes ever smaller; it decreases from 0.1 to 10 -7 in the range of 
the plot. 

4 Closing remarks 

Using a method introduced by Batchelor [2], the speed of propagation of a semi-infinite 
cylindrical bubble into a vortex has been calculated without knowning the shape of the 
nose and found to be comparable to the maximum rotational speed in the vortex. 

If the pressure in the reservoir is indeed constant in time, then the main assumption 
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in the present treatment is that of a semi-infinite bubble. This allows neglect of axial 
derivatives at station 2. How valid is this for a finite bubble of length LI. Since the axial 
velocity must vanish at the axial center of the bubble, axial derivatives of velocity ~ U\>/L. 
On the other hand, radial derivatives of velocity ~ ft hence axial derivatives are small 
provided U^/L <C 1. Since the numerical results showed that U b = 0(1), the requirement 
is that L 1. 

Using the assumption of quasi-steadiness, we applied the constant pressure solution to 
bubbles with variable pressure. The validity of the assumption can be checked a posteriori 
for each case by calculating the characteristic acceleration (34). Depending on physical 
parameters and initial conditions, three distinct cases are possible. They are summarized 
in the abstract. 

Let us list the most salient of the many questions that remain unanswered: 

(1) What happens when one tries to inject a fixed pressure bubble at a point (^Apb,T^ 
in the parameter space beyond a limit point? In a similar vein, what happens when a 
quasi-steady variable pressure bubble arrives at a limit point in scenario Bl. 

(2) What is the status of the gaps in figure 3? 

(3) The present work assumes that one already has a cylindrical bubble. What conditions 
make possible the evolution from a spherical to a cylindrical bubble? What determines 
whether the solution will end up on a type A or type B1/B2 trajectory. 

(4) A uniform axial velocity in the vortex cores merely adds to the velocity of bubble 
propagation. It would be useful to implement the present procedure (numerically) for more 
general profiles perhaps in conjunction with an experiment that injects air into a vortex 
ring or into the vortex of an airfoil in a water tunnel or towing tank. 

(5) The stability of a cylindrical bubble surrounding by a purely circumferential flow has 
been studied by Pedley [6], for instance. The result is that the surrounding flow stabilizes 
the bubble. In the prese nt situation the flow at station 2 also includes an axial flow. What 
are the stability properl ies of the station 2 flows? 
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